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In this clarifying note we discuss the late decay of an Affleck-Dine condensate by providing a no-go 
theorem that attributes to conserved global charges which are identified by the net particle number 
in fields which are included in the flat direction(s). For a rotating condensate, this implies that: 
(1) the net baryon/lepton number density stored in the condensate is always conserved, and (2) the 
total particle number density in the condensate cannot decrease. This reiterates that, irrespective 
of possible non-perturbative particle production due to D-terms in a multiple flat direction case, 
the prime decay mode of an Affleck-Dine condensate will be perturbative as originally envisaged. 
As a result, cosmological consequences of flat directions such as delayed thermalization as a novel 
solution to the gravitino overproduction problem will remain virtually intact. 



I. INTRODUCTION 

The scalar potential of the Minimal Supersymmetric 
Standard Model (MSSM) has many flat directions 
These directions are classified by gauge-invariant mono- 
mials of the theory, and most of them carry baryon 
and/or lepton number 0, The flat directions have 
many important consequences for the early universe cos- 
mology pj]. Most notably, there are two flat directions 
which can potentially act as the inflaton and can be 
tested at the LHC [1 (see also ^) \ 

Moreover it is well known that a baryon/lepton car- 
rying flat direction can generate the observed baryon 
asymmetry via Affleck-Dine mechanism ^1^. During in- 
flation a condensate is formed along the flat direction. 
After inflation, the condensate starts rotating once the 
Hubble rate drops below its mass. This results in a 
baryon/lepton asymmetry which will be transferred to 
fermions upon the decay of the condensate. The Vac- 
uum Expectation Value (VEV) of the condensate induces 
large masses to the flelds which are coupled to it. The 
decay to these fields will be only possible when the con- 
densate VEV has been redshifted to sufficiently small val- 
ues. This will result in a late perturbative decay of the 
fiat direction condensate [lo!| . A late decay of an Affleck- 
Dine condensate has another important consequence in a 
supersymmetric universe, na mely late thermalization of 
the inflaton decay products [III IT3 | . The flat direction 
VEV breaks the Standard Model (SM) gauge symmetry, 
thus inducing large masses to gauge (and gaugino) fields 
via the Higgs mechanism. This will slow down thermal- 
ization by suppressing dominant reactions which estab- 
lish kinetic and chemical equilibrium among the inflaton 
decay products ^. A delayed thermalization results in 



^ Other important applications of flat directions include curvaton 
mechasnism [^, inhomogeneous rehetaing and density perturba- 
tions [3, magnetic field generation @] , and non-thermal dark 
matter [^. 

^ Finite temperature effects on MSSM flat directions have been 



a reheating temperature much lower than what usually 
thought. This naturally solves the outstanding problem 
of thermal gravitino overproduction in supersymmetric 
models iUl] ^. 

The aim of the note is to underline the crucial im- 
portance of conserved global charges, which was first ob- 
served in seminal papers by Affleck and Dine [10 ] , and by 
Dine, Randall and Thomas Charges identified by the 
net particle number in fields which are included in a fiat 
direction, most notably baryon and lepton number, are 
preserved by the D-terms For a (maximally) rotating 
condensate, this implies that possible non-perturbative 
effects cannot change the baryon/lepton number density 
stored in the condensate, and will not decrease the total 
number density of quanta in the condensate. As we will 
briefly mention, under general circumstances, this also 
holds when _F-terms are taken into account. The decay of 
a rotating condensate into other degrees of freedom hap- 
pens through the F-term couplings. As discussed in the 
original work of Affleck and Dine , this decay occurs 
late and is perturbative ^. This guarantees that cosmo- 
logical consequences of an Affleck-Dine condensate, such 
as late thermalization, will proceed naturally. 

Similar conclusion arises for multiple flat directions 
represented by a gauge-invariant polynomial (for a de- 
tailed discussion, see Ref. [H, H^), as it is just a mani- 
festation of the conservation of global charges carried by 
a rotating condensate. 



discussed in [T^. Il4|. 

^ Thermalization in the presence of supersymmetric flat directions 
was first considered in [Tsl l (fo r some of the works in the non- 
supersymmetric case, see The effect of Q-balls formed 

from fragmentation of flat direction oscillations on reheating is 
discussed in ITtI . 

* The A-term does not preserve these charges. However it becomes 
irrelevant after the very first oscillations, as it is redshifted away 
rapidly 

^ It is known that the _F-terms cannot lead to a non-perturbative 
decay of a rotating condensate [24l | . 
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II. ROTATING FLAT DIRECTIONS 

A. Brief introduction to flat directions 

The scalar potential of the MSSM has a large number 
of flat directions. The D-term and F-term contributions 
to the potential identically vanish along these directions. 
The D-Aat directions are categorized by gauge-invariant 
combinations of the MSSM (super)fields $i. The D- 
flatness requires that ^ 



= 

i 



(1) 

(2) 



The rotation builds up already within the first Hubble 
time The ^-term is quickly redshifted compared 

with the mass term as it is higher in order (i.e. is a 
non-renormalizable term). It becomes negligible after a 
few Hubble times, thus leading to a freely rotating con- 
densate. The trajectory of the motion in the (j) plane is 
an ellipse, and its eccentricity is ~ 1 as the A-term is 
initially as large as the mass term 0]. For our purpose, 
it can be approximated with a circular trajectory: 



(pR = if cos{'mt) = if sin(mi). 



(4) 



where ip is redshifted cx a ^ due to Hubble expansion (a 
being the scale factor of the universe). 



T'^ are generators of the SU{3)c and SU (2) l symmetries, 
and Qi are the charges of under U{1)y- 

A subset of Z?-flat directions are also i^-flat, in a sense 
that the superpotential makes no contribution to the 
potential at the renormalizable level. However the flat 
directions are lifted by supersymmetry breaking terms, 
and non-renormalizable superpotential terms induced by 
physics beyond the standard model. Hence the potential 
along a flat direction, denoted by (j), generically follows 



M2("-3) 



AX 



h.c. 



(3) 



Here m, A ^ 0{TeV) are the soft supersymmetry break- 
ing mass and A-teim respectively. M is a high scale 
where new physics appears (like Alp or Mqut), « > 4, 
and A ~ C(l) typically. 

The flat direction fleld (p acquires a large VEV during 
inflation as a result of the accumulation of quantum fluc- 
tuations in that epoch. This leads to the formation of 
a condensate along the flat direction fi\ ^. After infla- 
tion the VEV of the condensate slowly rolls down. This 
continues until the time when the Hubble expansion rate 
is ~ TO, see Eq. The VEV of the condensate at this 
time is ipo ~ (mM"^^ )i/"-2 ^ ^ 2], and all three terms 
in Eq. ([3]) are comparable in size. 

Of particular importance is the A-term which, by ex- 
erting a torque, results in the rotation of the condensate. 



Here we use the same symbol for a superfield and its scalar com- 
ponent. 

The supergravity corrections to the flat direction potential are 
given by cH?^f\ip\'^ (^^inf being the Hubble expansion rate 
during inflation). IfO<c~C'(l), the flat direction will settle 
at the origin during inflation and will play no dynamical role af- 
terwards. However, if < c ■C 1, or for c < 0, the flat direction 
acquires a large VEV during inflation and plays an important 
role in the post-inflationary universe. We will be mainly con- 
cerned with the latter case. Note that if there is a positive 0(1) 
Hubble correction to the flat direction mass, then it would also 
be reflected in the inflationary potential which would spoil the 
success of inflation [J. 



B. Physical degrees of freedom 

Let us consider the simplest flat direction represented 
by the HuL gauge-invariant combination. Here is the 
Higgs doublet which gives mass to the up-type quarks, 
and L is a left-handed (LH) lepton doublet ^. After im- 
posing the Z?-flatness condition in Eqs. (|ll2p . one can 
always go to a basis where the complex scalar field (the 
superscripts denote the weak isospin components of the 
doublets and R, I denote the real and imaginary parts 
of a scalar field respectively) 



V2 



(5) 



represents a flat direction. The Vacuum Expectation 
Value (VEV) of 0, denoted by ip, breaks the SU{2)l x 
U{1)y down to ?7(l)om (in exactly the same fashion as 
in the electroweak symmetry breaking) . The three gauge 
bosons of the broken subgroup then obtain masses ~ gip 
{g denotes a general gauge coupling). 
After making the following deflnitions: 



Xi = 

X2 = 
X3 



V2 



(6) 



we find the instantaneous mass eigenstates 



X2it) 



cos{mt)xi.R + sm{mt)x 



i.i 



V2 

cos(mt)x2 + sin(mt)x3 



(7) 



V2 



The situation will be similar for the H^H^ flat direction. 
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which acquire masses equal to those of the gauge bosons 
through the D— term part of the scalar potential has 
both real and imaginary parts). Note that 

cos(mt)Y^ — sin(mt)Y^„ 

xm = ' (8) 



X4(i) = 



V2 

cos(mt)x3 — sm{mt)x2 



are the three Goldstone bosons (again Xa has both real 
and imaginary parts), which are eaten up by the massive 
gauge fields via the Higgs mechanism. Therefore, out of 
the 8 real degrees of freedom in the two doublets, there 
are only two physical light fields: (jjR and (/)/, i.e. the real 
and imaginary parts of the flat direction field ^ . 

A rotating fiat direction, see Eq. does not cross 
the origin. Hence, starting with a large VEV such that 
gip 3> TO, the hierarchy between the mass eigenvalues of 
the heavy and the light degrees of freedom is preserved 
at all times. However, rotation results in time variation 
in the mass eigenstates of the fields Eq. (I7I8|) . 

The heavy fields, despite having time- varying mass 
eigenstates ([7]), evolve adiabatically at all times since 
gif ^ TO, and hence will not experience any non- 
perturbative effects. In fact, they get decoupled and be- 
come dynamically irrelevant. 

If there are light fields with a mass < to, time vari- 
ation will become non-adiabatic and can lead to non- 
perturbative production of their quanta [13, [U, [l^, H^l . 
In the case at hand the only physical light degrees of free- 
dom, i.e. degrees of freedom on the D-flat subspace after 
removing the Goldstone modes ^ , are the real and imag- 
inary parts of the fiat direction field ([5]) . The mass eigen- 
states and eigenvalues of the flat direction are both con- 
stant in time, therefore, there will be no non-perturbative 
effects whatsoever (for detailed discussion, see Ref. [11]). 

In general, the number of physical degrees of freedom 
on the £)-flat subspace is given by [l^ 



Might = 



total 



(2 X A^brokcn)- 



(9) 



.^otai is the total number of degrees of freedom in the 
gauge-invariant combination which represents the flat di- 
rection(s), and A^'broken is the number of spontaneously 
broken symmetries. The factor 2 counts for the Gold- 
stone bosons which are eaten by the gauge flelds plus the 
superheavy degrees of freedom which are decoupled. 

In the case of H^L flat direction, Eq. ^ reads: 
Niigu = (2 X 2 X 2) — (2 X 3) = 2, as exphcitly shown above. 



^ The masses induced by the flat direction VEV are supersym- 
metry conserving. One therefore finds the same mass spectrum 
in the fermionic sector as in above (for details, see [isj'l. Note 
however that scalars also acquire soft supersymmetry breaking 
masses ~ C'(TeV), while fermions do not. This implies that the 
fermionic partner of the flat direction field (f) will remain exactly 
massless. 



This is the typical tendency for a single flat direction in 
MSSM [is], i.e. when the flat direction is represented by 
a gauge-invariant monomial. Therefore particle produc- 
tion due to rotation can only be possible if one considers 
two or more flat directions, typically represented by a 
gauge-invariant polynomial [2^ . 

Nevertheless, we would like to emphasize that time- 
variation in the mass eigenstates of light fields, even if 
it happens, has practically no bearing for the decay of 
rotating flat direction(s). This is the topic which we will 
discuss in the next section. 



III. NO-GO THEOREM FOR ROTATING FLAT 
DIRECTIONS 

Let us consider MSSM flat direction(s) represented by 
a gauge-invariant combination of the fields ^i. The D- 
term part of the potential, see Eqs. (|ll2p . is invariant 
under arbitrary phase transformations of (the same is 
true for kinetic terms) ^": 

$«^e'"'$,. (10) 

The associated conserved charges ("•" denotes differenti- 
ation with respect to time) 



(11) 



represent the net particle number density, i.e. the differ- 
ence between the number density of particles and anti- 
particles, in <f>i ^^: 



(^particle ^anti — particle)^ 



(12) 



The total particle number density in <I>i, denoted by fii, 
is the sum of particle and anti-particle number densities: 



'^i — (^particle ^anti — particle) j ^ \ ■ 



(13) 



For maximally rotating fiat direction(s), the fields <I>i fol- 
low: 



^/2 



exp(i6'i), 



(14) 



where (pi = 0, and 6i = from the equations of motion 
Eq. pT|) then reads 



(15) 



^'^ Some combinations of these phases are associated with the U{1)y 
and U{1) subgroups from diagonal generators of SU{3)c and 
SU(2)i^. Hence in a background with non-zero flat direction 
VEV they correspond to Goldstone bosons which, in the unitary 
gauge, are completely removed from the spectrum. We will deal 
with this carefully in the case of explicit examples presented in 
the next section. 

The j4-term, see Eq. JSj, breaks these symmetries. However it 
virtually decouples within a few Hubble times after the rotation 
starts. In particular it will be irrelevant by the time possible 
non-perturbative effects become important. 

For a freely rotating scalar field with mass rrii we have di = m^. 
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On the other hand, hi is given by (the factor of 2 accounts 
for the kinetic term and the mass term of (f>i) 

= 2|^,rV.|' = l^»U?- (16) 

The inequahty in Eq. (|13p is therefore saturated, hi = 
I nil, for maximal rotation. This is not surprising as such 
a condensate consists of particles or anti-particles only. 
This leads to a following no-go theorem: 

Consider MSSM flat direction(s) represented by 
gauge-invariant combinations of fields $i. Possible 
non-perturbative particle production from time-variation 
in the mass eigenstates caused by the D-terms: 

(1) cannot change the net particle number density 
in denoted by Ui, and hence the total baryon/lepton 
number density stored in the condensate. 

(2) cannot decrease the total particle number den- 
sity in $i, denoted by hi, thus the total number density 
of quanta h = '^ihi in the condensate . As a direct 
consequence of the conservation of energy density, 
non-perturbative effects will not increase the average 
energy of quanta E'avc • 

We note that, so far as the D-terms are concerned, 
the theorem also applies to the subsequent evolution of 
the plasma formed after the phase of particle production. 
This implies that possible non-perturbative effects do not 
lead to the decay of a rotating condensate. They merely 
redistribute the energy which is initially stored in the 
condensate among the fields on the D-flat subspace 

Some comments are in order before closing this sec- 
tion. Particle production drains energy from the rotat- 
ing condensate. The question is how this energy transfer 
affects the flat direction(s) trajectory. The I?-terms, see 
Eqs. (jll2p . and kinetic terms are invariant under inter- 
changing the real and imaginary components of scalar 
fields $i which are included in the flat direction(s). For 
a circular motion the trajectory itself is invariant un- 
der such interchanges. This implies that possible non- 
perturbative particle production, which is governed by 
D-terms and kinetic terms, will not change the shape of 
a circular trajectory. Therefore all that happens in the 
case of maximal rotation, see Eq. (jl]) , is a decrease in the 
radius of the circle if. 

In reality, the condensate will not undergo maximal 
rotation and trajectory of its motion will be an ellipse: 

4>fi = ip cos(mt) — aif sin(mi), (17) 



rather than a circle where a < 1 is a positive number. 
However, for a generic initial condition, the rotation will 
be near maximal a ^ 0.3 f^. The difference between 
the number density of particles and antiparticles in the 
condensate is given by 

n = 2am<y9^ (18) 

while the total number density of quanta in the conden- 
sate follows 

h = {I + a^)mip^ . (19) 

Note that for a = 1 the trajectory is a circle and n = h 
as we discussed above. However, for a ~ 1 we have n ^ h 
(for a < 1 we always have h > n). The condensate in this 
case consists mainly of particles (or anti-particles), but 
it also contains a small mixture of anti-particles (or par- 
ticles). Therefore, in agreement with the conservation 
of net particle number density, h can in principle de- 
crease by a factor of r = (1 + a'^)/2a, such that the small 
mixture of (anti-)particles will vanish and consequently 
h = n. For a ~ 0.3, we have r ^ 2. The possible decrease 
in h will therefore be of 0(1). As we will discuss later 
on, the situation is similar to that for a maximal rota- 
tion with regard to the final decay of the flat direction(s) 
energy density and late thermalization of the universe. 

Finally, we note that the i^-terms, due to quark 
and lepton mixing, preserve Ui (i.e. the total 
baryon/lepton number density) instead of each n^. How- 
ever, I J2i "-il ^ J2i I'^ili except in some cases where J2i "^i 
is much smaller than the the individual n^. This requires 
special initial conditions for which the baryon/lepton 
number density stored in individual fields are large, but 
come with opposite signs in such a way that that they 
conspire to make the total baryon/lepton number den- 
sity which is stored in the condensate much smaller 
Hence, under general circumstances, the no-go theorem 
holds when all interactions in the MSSM Lagrangian are 
taken into account. 



IV. SOME EXAMPLES 

To elucidate the no-go theorem, we consider three rep- 
resentative examples of MSSM flat directions. Namely, 
single flat directions consisting of two and three fields, 
and multiple flat directions. 



Here we mean the comoving quantities as the Hubble expansion 
inevitably redshifts any physical number density. 
As we will explain later, this is in sharp contrast to non- 
perturbative particle production from an oscillating condensate, 
also called preheating, studied in the context of infiaton de- 
cay 



The only case where this can happen naturally is for the HuHd 
flat direction which carries B = L = 0. However, this is a single 
flat direction for which there is no time variation in the light 
physical degrees of freedom [Tsl l . Therefore there can be no non- 
perturbative particle production in this case at the first place. 
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HuL direction 



The £)-terms associated with SU{2)l and ^7(1)^, see 
Eqs. ()ll2p . are invariant under two C/(l) symmetries: 



(20) 



and the corresponding charges 



ni = iH*Hu 



h.c. 



n2 



iiU + h.c, (21) 



are conserved. 

In a background of rotating flat direction, transforma- 
tions generated by non-diagonal generators of SU{2)l 
can be used to situate the VEVs along and (su- 
perscripts denote the weak isospin components), which 
we denote by 0i and 02 respectively: 

01 = I exp(i0i) , '/'2 = I exp(z02) . (22) 

The phase difference 6*2 — 0i is a Goldstone boson which 
can be removed through a U{\)y transformation (for 
identification of Goldstone modes, see Section IX). This, 
as shown before, leaves us with only two light degrees of 
freedom 



t 
2 



exp(i0 
Eq. (PT|) then results in 



62 = I exp(i0). (23) 



"1 



n2 = ip . 



(24) 



For a rotating flat direction ip — 0. Then from the equa- 
tions of motion we find 9^ = {rri^j + m\)/2, where mn 
and rriL are the masses of H and L respectively. Note 
that n2 is the lepton number density stored in the con- 
densate. 

The total particle number density in and L (de- 
noted by fii and ri2 respectively) follow from Eq. ([T| 



ni 
n2 



ip = \ni\ 
2 



\n2\ 



(25) 



and the corresponding charges 

ni = iu*Ui-\-\i.c. n2 = id* dj +h.c. 77.3 = id^dfc -|- h.c, 

(27) 

are conserved. 

In a rotating flat direction background, transforma- 
tions generated by non-diagonal generators of SU (3)c can 
be used to situate the VEVs along uj, (which we 

denote by 0i, 02, 03 respectively), where superscripts 
denote the color indices: 

01 = ^ exp(z6'i) , 



b2 ^ -j= exp(z6'2) , 



exp(i6'3). 



V6 



(28) 



The phase differences {29i — 6*2 — ^3) and 0i — 62 are 
Goldstone modes which can be removed through trans- 
formations generated by diagonal generators of SU{3)c ^'^ 
(for identification of Goldstone bosons, see Section XI). 
After the removal of Goldstone bosons, only two light 
degrees of freedom remain 



V6 



exp(i6 



exp(ii6 

ip 

71 exp(i6 



(29) 



For a rotating flat direction we have 1^ = 0. Then from 
the equations of motion we find 0^ = {m^. + rn^, + 



)/3. Eq. ^ resuhs in 



ni 



n2 = ris = 6 (p 



(30) 



Note that n = ni + n2 + is three times the baryon 
number density stored in the rotating condensate {u and 
d have baryon number —1/3). 

The total particle number density in Ui, dj, dk (de- 
noted by fii, 77,2, ns respectively) follow from Eq. (|16p : 



B. udd and LLe directions 

The situation for udd and LLe fiat directions is quite 
similar. We therefore concentrate on the udd case. 
The £)-terms associated with SU{3)c and U{1)y, see 
Eqs. (jll2p . are invariant under three U{1) symmetries 
(subscripts are the family indices): 



e'"'dk, 



(26) 



Or, equivalently, the diagonal generator of SU{2)l. 



hi 


= 1^1 = 


\ni 


n2 


= 1^1 = 


|"-2 


ns 


= 1^1 = 


|"-3 



(31) 



The action of U{1)y is the same as that of the {—1,-1, -(-2) 
diagonal generator of SU{3)c- 

This is different from the toy example presented in Ref. |22( . 
which considers a flat direction consisting of three fields charged 
under a single U{1) gauge symmetry. In the case of MSSM, 
there are enough symmetries to rotate away all phase differences 
among the fields, and hence only the overall phase remains as a 
physical degree of freedom. 
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C. HuLi multiple flat directions 

Now we consider multiple flat directions represented 
by the H^Li polynomial where all three Li dou- 

blets have a non-zero VEV. This case was first consid- 
ered in Ref. Q. The D-terms associated with SU{2)l 
and U{1)y, see Eqs. (|1I2I) . are invariant under four C/(l) 
symmetries: 



Li 



o'"3 r 
2 ^3 



(32) 



L2 ^ e*"^L2 

and the corresponding charges 
ni — iL\Li + h.c. 712 = iL\L2 

are conserved. 

Similar to the case of the single flat direction, the 
VEVs can be situated along the L\, L^, L^, compo- 
nents which we denote by 4>i, (j)2, 4>3^ 't>i respectively 



h.c. 713 ~ iL'^L^ + h.c. , 
h.c. , (33) 



19. 



2 

— exp{i02) 
2 

^ cxp(i6l4) , 



exp(i0i 



exp(«6l3 



(34) 



where Lp^ ~ ipl + + (fi"^ is imposed by the D-flatness 
condition, see Eqs. (|ll2p . The phase {(p9 — (pi9i — (/32^2 — 
tpzOs) is a Goldstonc mode which can be removed by a 
U{1)y transformation 2° (for identification of Goldstonc 
bosons, see Section X). After its removal we can recast 
Eq. ([34]) in the following form 



exp(i0i) 
— exp{t02) 

2 
f_ 
2 



exp(i6l3) , 



exp 



(35) 



Eq. ([33]) now results in 



ni = 






n2 = 






ns = 


^3 ^l, 




714 = 


{ipiOi + ip202 + ^3^3) V- 


(36) 



For maximal rotation (pi, ip2, (fa are constant, and 9i = 
02 ~ S3 ^ from the equations of motion 

Note that 71 = rii -1-7^2-1-77-3 is the Icpton number stored 
in the condensate. The total particle number density in 
Li, L2, L3, Hu (denoted by 77.1, 77.2, 77,3, 774 respectively) 
follow from Eq. (fT6|) : 



771 = 




= |77l 


772 = 


\02\'A 


= |772 


773 = 


\03\^\ 


= |773 


77.4 — 


- 


f </52^: 



(37) 



V. DIFFERENCES BETWEEN ROTATING AND 
OSCILLATING CONDENSATES 

Let us consider an oscillating condensate for which the 
trajectory of motion is a line instead of a circle: 



— cos{mt) . 



= 0, 



(38) 



where g(p 3> (g is a typical gauge coupling). In this 
case the mass eigenstate of the x fields which are coupled 
to 4> through the ZJ-terms, see Eq. ([B]), are constant in 
time but the mass eigenvalues oscillate. The time vari- 
ation becomes non-adiabatic as ~ 0. As a result, x 
quanta are created within short intervals each time that 
(j) crosses the origin [26] . This leads to an explosive stage 
of particle production, also called preheating, which even- 
tually results in a plasma of x and (j) quanta with typical 
energy 



(gLpm) 



1/2 



> m. 



(39) 



This implies an increase in the average energy of quanta, 
and hence a decrease in the number density of quanta, 
as compared to the original condensate. If the universe 
were to fully thermalize after preheating, we would have 
Save ^ T ^ {mipY^'^. Preheating is therefore a step 
toward full thermal equilibrium as it partially increases 
i?avo toward irs equilibrium value 

This is in sharp contrast to the situation for a rotating 
condensate. There, as we argued, possible particle pro- 
duction cannot decrease the number density of quanta. 
The marked difference between the two cases can be un- 
derstood from the trajectory of motion (i.e. circular for 
rotation versus linear for oscillation) . An oscillating con- 
densate 6 can be written as 



- exp (76 



2 exp( 



(40) 



The situation is actually more subtle than that for the HuL single 
direction, see Section IX. 

Or, equivalently, a transformation generated by the diagonal gen- 
erator of SU{2)l. 

Note that if any two of ipi, ip2, are zero, the situation will 
be reduced to that for the HuL single flat direction in Eq. II23I I. 



The situation is more complicated than the previous examples 
due to having more than one physical phase. In particular, ipi 
can rotate with different velocities, where the velocity of rotation 
is given by 8i. 

In reality, it takes a much longer time to establish full thermal 
equilibrium after preheating [2711 . 
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and the conserved charge associated with the global U{1) 
(corresponding to phase 9) is given by 

71 = i0*0 + h.c. = 0. (41) 

This is not surprising since an oscillation is the superpo- 
sition of two rotations in opposite directions, which carry 
exactly the same number of particles and anti-particles 
respectively. Therefore the net particle number density 
stored in an oscillating condensate is zero. 

Now consider non-perturbative particle production 
from an oscillating condensate. One can think of this 
process as a series of annihilations among N particles 
and N anti-particles in the condensate, > 1, into an 
energetic particle-anti-particle pair. This is totally com- 
patible with conservation of charge, see Eq. (|4T|) : n = 
after preheating as well as in the condensate. 

On the other hand, a (maximally) rotating con- 
densate consists of particles or anti-particles only, see 
Eqs. (P51[5n[57|) . Conservation of the net particle number 
density then implies that N ^ 2 annihilations {N > 2) 
are forbidden: annihilation of particle (or anti-particle) 
quanta cannot happen without violating the net particle 
number density. Therefore the total number density of 
quanta will not decrease, and the average energy will not 
mcrease . 



VI. DECAY OF A ROTATING CONDENSATE 

As we have discussed, any possible non-perturbative 
particle production will result in a plasma which is at 
least as dense as the initial condensate. All that can hap- 
pen is a redistribution of the energy density in the con- 
densate among the fields on the D-flat subspace. These 
fields have masses comparable to the flat direction mass 
m, as they all arise from supersymmetry breaking. Then, 
since the average energy is -Eavc < fn, the resulting 
plasma essentially consists of non-relativistic quanta. Its 
energy density p — ri-Eavo is therefore redshifted cx 
(a is the scale factor of the universe). 

The question is when this plasma will decay to other 
MSSM fields, in particular fermions, and thermalize. The 
plasma induces a large mass rrics to the scalars which are 
not on the Z?-flat subspace and their fermionic partners 
through the F-terms. In the Hartree approximation the 
effective mass is given by (for example, see [28j ) 

mls-h'l^, (42) 



where h denotes a Yukawa coupling ■^^ . The one-particle 
decay is kinematically forbidden as long as > m. 
Note that higher order processes such as iV ^ 2 anni- 
hilations (TV > 2) cannot happen due to conservation of 
global charges (i.e. baryon and lepton number) in the 
plasma. Since n and i?avc are respectively > and < than 
their corresponding values in the initial condensate, thus 
TOoff will always be larger than the induced mass by the 
condensate VEV, which is given by hip. 

Further note that rUeS is redshifted (x a^'^/^, where a cx 
//-2/3 j-^-i/2-| jj-^ matter (radiation) dominated epoch. 
The decay of energy density p (initially stored in the 
condensate) happens only when m^s has been redshifted 
below m, at which time the Hubble expansion rate is 
given by {ipo is the initial VEV of the condensate) (Tll.[l^ 

Hdcc ~ ni I (matter domination) , 

\hipoJ 

f m V'^ 

Hdcc ^ nil ) (radiation domination) .(43) 

\h(poJ 

Hence for a large ipo the decay time scale is suffi- 
ciently large compared to the time scale for possible non- 
perturbative particle production, which is ~ m~^. This 
decay happens perturbatively as discussed in 0, • 

As we discussed earlier, in reality the condensate has 
an elliptic trajectory whose eccentricity is ~ 0.3. This, 
see Eqs. (|18|19|) . implies that h can at most decrease (and 
hence i?avc increase) by a factor of 2 compared with their 
corresponding values in the initial condensate. Therefore 
TOeff, see Eq. (|42|) . may be smaller by a factor of 2 in 
the case of a realistic elliptic trajectory. According to 
Eq. this will result in an i/dcc which is larger by a 

similar factor factor, thus a slightly earlier perturbative 
decay. Nevertheless, for ipQ^ m, the final decay happens 
much later than the initial phase of non-perturbative par- 
ticle production. 

This reiterates the main point of this note: a phase of 
non-perturbative particle production due to rotation, al- 
though possible, cannot lead to the decay of fiat direction 
energy density. It will merely result in a redistribution of 
energy density on the D-fiat subspace. The final decay 
(to other fields) will happen late, and will be perturba- 
tive, as originally envisaged 0, HO]. 

VII. COSMOLOGICAL CONSEQUENCES 

For practical purposes, the resulting plasma will be- 
have the same as the initial condensate. In this section we 
discuss some of the important cosmological consequences 
for an Affleck-Dine condensate which gives rise to delayed 
thermalization and a solution to the gravitino problem. 



Note that an increase in the total particle number density, 
through creation of an equal number of particles and anti- 
particles will be in agreement with the conservation of the net 
particle number density. In this case the resulting plasma will 
be even denser than the condensate. 



In the case of thermal equilibrium (and zero chemical potential) 
we have n ~ and Save ~ T, where T is the temperature, 
yielding the familiar result ~ h?T^ . 
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A. Delayed thermalization 

The condensate VEV, denoted by spontaneously 
breaks the SM gauge symmetry and induces a large mass 
for the gauge/gaugino fields rrics ~ 9^ via the Higgs 
mechanism. Such a large mass suppresses gauge inter- 
actions which play the main role in establishing thermal 
equilibrium among inflaton decay products [11]. For a 
rotating condensate (p changes only due to the Hubble 
redshift. The gauge interactions will therefore remain in- 
effective for a long time until ip has been redshifted to a 
sufficiently small value. It is only at this time that full 
thermal equilibrium can be established [HI, ■ 

Now consider the plasma consisting of quanta of the 
fields <i>i on the D-flat subspace. The SM gauge symme- 
try is broken in the presence of this plasma as well. This 
results in an induced mass mcS for the gauge fields (and 
gauginos) which, as mentioned earlier, is given by 



mis 
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(44) 



Since n and -Eavo are respectively > and < than the cor- 
responding values in the initial condensate, it turns out 
that TTieff > 9ip. This implies that the gauge interactions 
will be (at least) as suppressed as that in the presence 
of a condensate. Hence, considering that the plasma de- 
cays like the initial condensate, thermalization will also 
be delayed similarly, for details see [ll|, . 

Again we note that for a realistic elliptic trajectory 
mefif may be smaller by a factor of 2. However, for ipo 3> 
m, universe thermalization will still be considerably de- 
layed relative to an initial phase of non-perturbative par- 
ticle production from the rotating flat direction(s). 



B. Thermal generation of gravitinos 

Late thermalization has an important consequence for 
thermal production of gravitinos [2^^^. 

First, delayed thermalization leads to a considerably 
low reheat temperature given by the expression [lli] 



(FthrA/p)'^' , 



(45) 



1 /"^ 

instead of the usual expression Tr ~ (FdMp) Here 
Fthr is the rate for thermalization of the inflaton decay 
products and Fd is the inflaton decay rate. Suppression 
of the interactions that lead to establishment of thermal 
equilibrium, due to the VEV of flat direction(s), implies 
that Fthr ^ Fd [llj, and hence a much lower Tr than 
usually found. 



Non-thermal gravitino production at early stages of inflaton os- 
cillations [30)1 is not a major issue as discussed in [3l| . 



In particular, one can naturally obtain Tr ^ 10^ GeV 
which, in the case of a weak scale supersymmetry, is re- 
quired in order not to distort predictions of the Big Bang 
Nucleosynthesis by the decay of gravitinos [29]. 

Moreover, before thermalization of the inflaton decay 
products, scattering processes which lead to gravitino 
production make a negligible contribution (for details, 
see [ll|). These two effects address the long standing 
gravitino problem in a natural way within supersymme- 
try without invoking any ad-hoc mecahnism. 



VIII. CONCLUSION 

The important message of this note is that possible 
non-perturbative effects stemmed from the D-terms have 
no bearing for the decay of energy density in rotating flat 
direction(s). This is due to the conservation of global 
charges associated with the net particle number density 
in flclds which are included in the flat direction(s), most 
notably the baryon/lcpton number density [lo|. For a 
rotating condensate, this ensures that the total number 
density of quanta will not decrease and, consequently, the 
average energy of quanta will not increase. 

Thus, in sharp contrast to an oscillating condensate 
(as in the case of inflaton decay via preheating) , all that 
can happen is a mere redistribution of the condensate 
energy among the fields on the D-fiat subspace. The 
actual decay into other fields happens pcrturbatively as 
originally envisaged by Affleck and Dine [10]. This en- 
sures the success of cosmological consequences such as 
delayed thermalization as a novel solution to the grav- 
itino problem [ll|, [l2| . 
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X. IDENTIFICATION OF GOLDSTONE 
BOSONS 

Here we quickly comment on identification of the Gold- 
stone bosons and their removal from the spectrum. For 
simplicity, we consider the case with a single U{1) gauge 
symmetry, but generalization to non-abelian symmetries 
is straightforward. 

Consider n scalar fields 4'i with respective charges qi 
(1 < i < n) under the U{1) symmetry. Covariant deriva- 
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tives of the scalar fields are 

n 

^ {d^ + tA^)cb* {d^ - (46) 
1=1 

The scalar fields can be written in terms of radial and 
angular components (denoted by (p and 9 respectively): 

= ^ exp(i0,)- (47) 

Expanding the fields around a background where ipi is 
constant, as happens for rotating flat direction(s), we 
then have 

= i{d^ei) exp(z0,)- (48) 

It can be seen that the combination X]"=i [fili^^i) can 
be eliminated from Eq. (|46p by performing the following 
gauge transformation: 

n 

^ A^-Y,<l^i.^^,e,). (49) 

1=1 

Therefore it is not a true physical degree of freedom. 
This particular combination is nothing but the Goldstone 
boson from spontaneous breaking of U{1) symmetry by 
non-zero values of Lpi. 

XI. Y.i HuLi MULTIPLE FLAT DIRECTIONS 

Consider a general VEV configuration of and Li 
(1 < i < 3) that satisfies the Z?-flatness condition in 



Eqs. pi2p . We can always use non-diagonal generators 
of SU{2)l to rotate iJ„ to a basis where {Hi) — (su- 
perscripts denote the weak isospin component). 

In the case of H^L single flat direction (where two 
of the Li have zero VEV), D-flatness under the non- 
diagonal generators directly implies (L^) = in this ba- 
sis. However, for multiple flat directions, it is not so 
obvious that (Lf) = (L^) = (Lf) = in the basis where 
(ffi) = 0. 

Hence let us consider a general configuration where 
both isospin components of Li have a non-zero VEV. 
The vanishing of the D-term associated with the diagonal 
generator of SU{2)l then implies 

3 3 

-|(i/^)P+^|(ii)p-^|(L,^)P=0, (50) 

i=l i=l 

while vanishing of the D-term associated with U{1)y re- 
quires that (note that Hu and Li have opposite hyper- 
charge quantum numbers) 

K^^5P-Ei(il>p-Ei(i?>p = o. (51) 

1=1 1=1 
It is readily seen that the third term on the LH side of 

Eqs. (|50l5ip must vanish, thus {Lfj = (L|) = {Lfj = 0. 

Also: 

i=l 
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